Asymptotic Dipole Radiation Formulas 

By W, HOWARD WISE 

THE analysis of the radiation from dipoles as given by Sommerfeld 
and by von Hoerschelmann is deficient in one respect: it does 
not give the true l asymptotic expressions for the radiation leaving at 
a considerable angle from the horizontal. The correct asymptotic 
formulas have already been easily supplied by an appeal to the 
Reciprocal Theorem; lately M. J. O. Strutt 2 has got them directly 
from the boundary conditions and H. Weyl 3 has derived the correct 
asymptotic formula for a vertical dipole at the surface of the earth 
by a method quite different from Sommerfeld's. In the present 
paper it is shown how they can be got by merely improving the rigor 
of Sommerfeld's analysis. 

The present analysis begins with the formulas of von Hoerschelmann 
for the wave potentials of vertical and horizontal dipoles at a finite 
distance above the surface of the earth and generally follows Sommer- 
feld. The derivation of an asymptotic approximation for the wave 
potential of a vertical dipole is considerably different from Sommer- 
feld's and results in the simpler and more precise formulas deduced 
from the reciprocal theorem. 

Most of the analysis is somewhat simplified by taking the permea- 
bility of the earth to be unity. 

The notation used is chiefly that of Bateman. 4 

t = variable of integration, throughout the paper. 

k 2 = €/ioi 2 + 4:ir<r/iiu 

I = Vt 8 - ki 2 , m = Vr 2 - k 2 2 . 

The subscripts 1 and 2 refer to air and ground respectively. 
Ri, i?2, a, p, f, 10, x, y and z are adequately defined by Fig. 1. 

cos X = x/R, cos By = y/R, cos 0, = z/R, R 2 = x 2 + y* + z 2 . 
The wave potential of a horizontal dipole is 5 > 4 

1 See paragraph following equation (8). 

2 M. J. O. Strutt, Ann. d. Phys., Bd. 1, p. 721, 1929. 

3 H. Weyl, Ann. d. Phys., Bd. 60, p. 481, 1919. 

4 "Electrical and Optical Wave Motion," pp. 73-75. 

6 H. v. Hoerschelmann, Jahrb. d. draht. Teleg., Bd. 5, 1912, pp. 14-188. 
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When k 2 = °° , V = 2i? 2 _1 exp ikiRa and £/ is zero. When 
k 2 = ki, V = lR 2 ~ y exp ikiRa and £ Z A just cancels the field of the 
image dipole in £ M . 

The wave potential of a vertical dipole is 5 - 4 

n --=V-V+ 7 - (4) 

7 is the function n analysed by Sommerfeld in Riemann-Webers 
Differentialgleichungen der Physik. Sommerfeld transforms the inte- 
gration from to co into an integration from — °o to + °° by re- 
placing the Bessel function by its equivalent Hankel functions and 
then wraps the real axis path of integration around the zero of k 2 2 l 
+ ki 2 m and the two branch cuts from ki and k 2 to + i <* . He thus 
gets V = P + Qi + Qi where P is the integral around the zero of 
k 2 -l + ki 2 m, Qi is the integral around the branch cut from k u Qi is 
the integral around the branch cut from k% and the function inte- 
grated is 

I v k fu 2 Ht{T P )e-» l TdT. (5) 

k 2 -I + kftn 

I and m are pure imaginaries on the branch cuts, which are rectangular 
hyperbolas, from ki and k 2 respectively. As one carries t counter- 
clockwise around the branch cut from ki to + i «> , / travels up the 
right hand side of the imaginary axis from — *oo + e to + i<x> + e. 
A similar statement holds for m. 

Multiplying numerator and denominator of equation (5) by k 2 2 l 
— k]?m and then taking out a factor (k 2 4 - &i 4 ) _l our integrand 
becomes 



k 2 2 M-J? - k x 2 - kSr 2 - k 2 2 



IIo l (T P )e- w ^-^rdT, 



W-kf (r -S)(r + S) 

where s = + kik 2 + "V&2 2 + ki 2 . Integrating around the pole at 
t = 5, we get 



P = 
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= if k 2 ^s 2 - fci 2 - kMs* - k 2 * = 

- - 2 ^{% g o , (^ <wt|,h . ( 6 ) 

ktff _ ki* - kHs 2 - k 2 2 = 2k 2 Us 2 - k^ = 2ik x k 2 s. 
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/ and m must have their real parts positive and so in taking the 
square roots of s 2 — ki 2 = — &i 4 -5- (k 2 2 + &i 2 ) and s 2 — k 2 2 = — k£ 
■*■ (&2 2 + &i 2 ) one halves the smallest angle with the positive real axis. 
If they both lie on the same side of the real axis P is zero. In order 
that they may lie on opposite sides of the real axis it is necessary that 

arg ki 4 < arg (ki 2 + k 2 2 ) < arg k 2 \ 

Writing k 2 = a + ifi and k 2 2 = x + iy this means that 

2a0 , a 2 + fl 2 „ 

y > ., 0-. -V+ .. «. &■ 

a- — p- a~ — p- 

The goal of the paper being asymptotic formulas for the sky waves 
of vertical and horizontal dipoles the ground wave, P, will hereafter 
be ignored. This is possible because at the high frequencies for which 
dipoles are useful the ground wave is very highly damped. 

Sommerfeld gets an asymptotic expression for Q\ by noting that if 
we are at a great distance from the source most of the value of the 
integral comes from that portion of the path of integration very close 
to k\. The solution he arrives at is 

_2(fl+$2' + fl'+ + •••)^ where - k f -|. 8 (7) 

Neglecting higher powers of \jR than the first, equation (7) sums 
up into 

r. 2k 2 2 cos 6 t e iklR ,_ N 

k 2 2 cos e, + kdk 2 ~ - k; 1 R 



But in getting equation (7) Sommerfeld has replaced Vt 2 — k 2 2 by 
V&i 2 — k 2 2 . This is a needless approximation which ruins the sym- 
metry, damages the utility and tends to hide the physical meaning of 
the final result. To get the true asymptotic formula for Qi it is 
necessary to confine the approximations to the purely operational 
variety, i.e. make no approximations of substitution before integrating 
but let the approximation reside wholly in the manner of integrating, 
as follows 

= - ^ fW(rp) — (C„ + C x l + C 2 l 2 + • • .)rdr 

*Annalen der Physik, Band 28, 1909, page 705. 
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where 

Co + c x i + cy 2 + • • • = k 2 2 + (k 2 n + &i 2 V/ 2 + k? - w). 

To the extent that most of the value of the integral comes from 
that portion of the path of integration very close to t 2 = ki 2 the 
expansion in powers of / is valid. Replacing each /by — (d/dw) 
we have then 

0.~-2 = (G-C 1 £+fti- + - , )*£- (9) 



dw\ dw dw 2 J R> 
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+ higher order terms 

2£ 2 2 7 e<*«* 



& 2 2 7 + *i 2 V7 2 + *i 2 - ^2 2 ^2 

/ d ikx sin 2 6, a 2 \ 2fe 2 2 7 e"' g ' 

+ ^ cos 0, dy + 2 dy2 J k ^ + k ^—-^-— R j 

+ higher order terms. (10) 

The second order term can be neglected if kiR» > > 1, say if 

i? > 20X. It will hereafter be supposed that this is the case. 

Multiplying numerator and denominator of the first order term by 
i and canceling a ki we have finally 
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q 2fe 2 2 cos d z g 'hft 

1 ~ k 2 2 cos 6 2 + &W&2 2 - £i 2 sin 2 e, ^2 

This Qi behaves as one would rightfully expect a true asymptotic 
formula to behave. It is 2(exp ikiR 2 )/R 2 at k 2 = °o and l(exp ik\.R 2 )l 
i?2 at & 2 = &i- 

In so far as & 2 is considerably larger than ki and the expansion in 
powers of / is valid Q 2 is negligible in comparison with Q\. n Perhaps 
the easiest way of seeing this is to note that equation (10) might just 
as well have been obtained directly from V instead of from Qi. 

Substituting equation (11) in equation (4) we get 

IV ~ ( 1 + *'' cos °' - k ^ ~ kl * SJke^*-™ ) e — (12) 

\ k 2 2 cos e. + &W& 2 2 - W sin 2 e a ) Ri ' 

whence 



£/= -*«( n/ + ^!- 2 — n/ 



as 2 



- ita sin 2 5,(1 + Rtf***"***') ^- , (13) 



where 

„ _ k 2 2 cos 9 Z — &1V&2 2 — ki 2 sin 2 d z 
1 ~ fc 2 2 cos B, + M&2 2 ~ £i 2 sin 2 6, 

Substituting equation (11) in equation (2) and adding E z m we get 

gikiB 

E z h ~ - iu cos 0* cos 0,(1 - _R ie "i2««Kw».)— — ■ (14) 

i?i is the coefficient of reflection for a plane wave polarized in the 
plane of incidence. 

The horizontal fields of a horizontal dipole are 

v--*[v + *.-£(£h*+£iv+£il')] 

and 

^^-^[nz + ^^n^ + Anz + ln/)] 

7 Riemann-Weber's " Differentialgleichungen der Physik," p. 556. 
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= - «'«(n x " + n x A ) + cot <pE y h ; (16) 

/. £/ = E x h sin <p - E y h cos <p ico sin ^(n x °° + H, A ). (17) 

Evidently the procedure which yielded the true asymptotic expres- 




Fig. 2 — Vertical dipole polar diagrams computed for X = 6 meters, 
6 = 9, ix = 1, and <r = lO - " 

sion for F will do the same for IT X A . The details are not interesting. 
The result is 

2^i cos 0, «**»** 



n. 



£i cos 6, + V&2 2 - ki 2 sin 2 Z -^2 



(18) 
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Substituting; equation (18) in equation (17) we get 

oikiR 



E 9 h ~ + iu sin <p(l - RtfiWnO*) -—, 
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(19) 



where 



R« = 



V^a 2 — kr sin 2 6 Z — k\ cos 0- 
■yjk-/ — A' r sin 2 6 Z + k\ cos d z 



Ri is the coefficient of reflection for a plane wave polarized perpen- 
dicular to the plane of incidence. 

Formulas (13), (14) and (19) are just what one woukTget by applying 

15° 




Fig. 3 — Polar diagrams of the horizontal fields of horizontal dipoles computed for 
X = 6 meters, e =* 9, ju = 1 and a = 10 -13 

the coefficient of reflection from an imperfect reflector to the reflected 
waves of the corresponding electrostatic formulas. 
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Returning to equation (15) 
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Usually one cares only for E 9 h at ^ = x/2 and £,* and £„* are of 
no particular interest. Figs. 2, 3 and 4 are polar diagrams in the 




Fig. 4— Polar diagrams of the vertical fields of horizontal dipoles computed for 
X = 6 meters, e = 9, /i = 1 and a = 10 13 

vertical plane of equations (13), (14) and (19). Assuming the con- 
ductivity of the air to be zero and the dielectric constant to be unity 

k-r = &r(e + i2c\a), ki = 2tt/X, 

where e is the dielectric constant of the ground referred to air as 
unity and a is the conductivity of the ground in electromagnetic units. 
The values of e and a used in computing the polar diagrams are 
generally supposed to be somewhere in the neighborhood of their 
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average values for earth but they vary so much with the locality that 
the diagrams can scarcely be regarded as giving more than a general 
idea as to what may be expected of the formulas. 

The number attached to each curve is the height of its dipole in 
quarter wave-lengths. 

Formulas (13), (14) and (19) are just what one would get by applying 
the Reciprocal Theorem to two dipoles, one near the earth and the 
other far away. The electric field acting on the one near the earth 
is composed of a direct field and a reflected field which is i?i or R 2 , 
as the case may be, times the direct field. 

When the depth to groundwater, bedrock, an orebody or any 
marked discontinuity in the electrical properties of the ground is 
known and is not too great the effect of this discontinuity on the polar 
diagram ought not to be ignored. The asymptotic formulas for any 
stratified ground are got by putting the coefficients of reflection for 
a plane wave reflected from the surface of that ground in place of the 
corresponding coefficients in equations (13), (14) and (19). For a 
number of rather obvious reasons it would usually be out of the 
question to deal with more than one plane of discontinuity; one is 
bad enough. The coefficients for a single plane of discontinuity at a 
depth A are 



and 



„ . k» 2 cos 6 Z — rjiki-yjk^ — k{' sin'- d z 

&i = , 

k-i 1 cos d, + J71&W&2 2 — k{ 2 sin 2 8 t 



where 



„ , _ V&o' 2 — k{ 2 sin 2 6 Z — 77261 cos d z 
V&2 2 ~ ki 1 sin- 6 Z — ?72&i cos 9 Z 

_ n« 1 + 61 __ Hz 1 + h 

Mi 1 — 01 Mi 1 — o 2 

. ._ kfux^kF- kr sin- Ol - fe 3 2 M2-\%2 2 - £i 2 sin 2 d z gA ^ fa8 _ tlt , ln8fl . 

\ — . - .... t.- ■ j 

&2V3V&3 2 — ki 2 sin 2 6, — &3V2V&2 2 — &1 2 sin 2 6 z 
and 

5 2 = M3-Vfe 2 2 ~ k? sin 2 "^: - wW^ k{> sin 2 ^ / , 12AN . ^,_ t , Jaln3tf: . 
ti3^k 2 2 — k{ 2 sin 2 6 Z -f /xo-yjkr — k{ 2 sin 2 d. 

If A is not large and kz is considerably different from k 2 then 771 
and 772 will differ considerably from unity. 



